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Abstract. 

We explore, by means of a large ensemble of SPH simulations, how the level of turbulence affects the collapse and fragmentation 
of a star-forming core. All our simulated cores have the same mass (5.4M Q ), the same initial density profile (chosen to fit 
observations of L1544), and the same barotropic equation of state, but we vary (a) the initial level of turbulence (as measured 
by the ratio of turbulent to gravitational energy, a tulb = t/ turb /|fi| = 0, 0.01, 0.025, 0.05, 0.10 and 0.25) and (b), for fixed ar tur b, 
the details of the initial turbulent velocity field (so as to obtain good statistics). 

A low level of turbulence (a tm ±, ~ 0.05) suffices to produce multiple systems, and as a- tlll -b is increased, the number of objects 
formed and the companion frequency both increase. The mass function is bimodal, with a flat low-mass segment representing 
single objects ejected from the core before they can accrete much, and a Gaussian high-mass segment representing objects 
which because they remain in the core grow by accretion and tend to pair up in multiple systems. 

The binary statistics reported for field G-dwarfs by Duquennoy & Mayor (1991) are only reproduced with a- tllr b ~ 0.05. For 
much lower values of a t urb (~ 0.025), insufficient binaries are formed. For higher values of a tm ±, (i 0.10), there is a significant 
sub-population of binaries with small semi-major axis and large mass-ratio (i.e. close binaries with components of comparable 
mass). This sub-population is not present in Duquennoy & Mayor's sample, although there is some evidence for it in the pre- 
Main Sequence population of Taurus analyzed by White & Ghez (2001). It arises because with larger a t urb, more low-mass 
objects are formed, and so there is more scope for the binaries remaining in the core to be hardened by ejecting these low- 
mass objects. Hard binaries thus formed then tend to grow towards comparable mass by competitive accretion of material with 
relatively high specific angular momentum. 
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1. Introduction 

Turbulence appears to play a crucial role in the structure and 
evolution of molecular clouds, in the formation of star-forming 
cores within molecular clouds, and in the collapse and frag- 
mentation of cores to form protostars. 

The main evidence for turbulence in molecular clouds 
comes from their apparently fractal substructure (e.g. 
Elmegreen & Falgarone 1996; Elmegreen 2002), and from the 
almost universal power-law scaling relations between size (L) 
line-width (Av) and mass (M) 1 and the almost universal power- 
law mass function 2 to which this substructure subscribes, over 
many orders of magnitude, from the largest molecular cloud 
complexes (~ 10 2 pc, ~ 1O 6 M ), down to the smallest re- 
solvable structures (~ 10~ 2 pc, ~ 10~ 4 M Q ) (e.g. Larson 1981; 
Myers 1983; Stutzki & Giisten 1990; Hobson 1992; Hobson et 
al. 1994; Williams et al. 1994; Elmegreen & Falgarone 1996; 
Kramer etal. 1996; Kramer et al. 1998; Heithausen et al. 1998). 



Send offprint requests to: Simon.Goodwin@astro.cf.ac.uk 

1 L oc M a with 0.35 <, a i 0.55 , and Av oc M 13 with 0.2 < /3 < 0.3 

2 dN/dM oc M~T with 1.5 < y < 2.0 



Until recently, it had been presumed that molecular clouds 
were long lived, being supported against collapse by their inter- 
nal turbulence, and this was advanced as the reason for the low 
overall efficiency of star formation.. However, it is now recog- 
nized that turbulence cannot support clouds for long, because 
- even with a frozen-in magnetic field - the turbulence dissi- 
pates on a dynamical timescale (Mac Low et al. 1998; Stone 
et al. 1998). Instead clouds are presumed to form and disperse 
on a dynamical timescale, without ever reaching equilibrium 
(Ballesteros-Paredes et al. 1999; Elmegreen 2000; Pringle et 
al. 2001; Hartmann et al. 2001). 

In this highly dynamical scenario, cores form wherever a 
sufficiently dense and coherent converging flow is created by 
the turbulent velocity field (Elmegreen 1997; Padoan et al. 
1997; Hartmann et al. 2001; Klessen & Burkert 2000, 2001; 
Klessen et al. 2000; Padoan & Nordlund 2002; Mac Low 
& Klessen 2004). Frequently these cores are not gravitation- 
ally bound, and therefore they disperse soon after they form. 
However, occasionally they are gravitationally bound, and in 
this case they are likely to proceed straight into gravitational 
collapse; these are the cores we identify as 'prestellar' . 
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This scenario is in contrast to the idea that prestellar cores 
are supported magnetically, and evolve quasistatically by am- 
bipolar diffusion, until they become magnetically supercritical 
and collapse (e.g. Basu & Mouschovias 1994, 1995a,b; Ciolek 
& Mouschovias 1993, 1994, 1995; Morton et al. 1994; Ciolek 
& Basu 2000). The main effects of the quasistatic ambipolar 
diffusion phase are (i) to give the core more time to lose angu- 
lar momentum by magnetic braking, (ii) to organize the mate- 
rial so that its subsequent collapse is rather well focussed, and 
(iii) to allow turbulence to decay. All three effects mean that 
such cores are less likely to form multiple systems. Since most 
stars are observed to be in multiple systems (e.g. Duquennoy & 
Mayor 1991; Fischer & Marcy 1992; White & Ghez 2001), and 
since there is no observational evidence for magnetically sub- 
critical cores (e.g. Crutcher 1999; Bourke et al. 2001; Crutcher 
at el. 2003), we shall assume that ambipolar diffusion does not 
play an important role in the evolution of prestellar cores. For 
simplicity, we ignore the magnetic field altogether. 

In the highly dynamic scenario the collapse of a prestellar 
core is far more likely to lead to fragmentation and the forma- 
tion of multiple systems (e.g. Whitworth et al. 1995; Turner 
et al. 1995; Whitworth et al. 1996; Klein et al. 2001, 2003; 
Bate et al. 2002a,b, 2003; Bonnell et al. 2003; Delgado-Donate 
2003, 2004; Goodwin et al. 2004a,b; Hennebelle et al. 2003, 
2004). This is because in the highly dynamic scenario prestel- 
lar cores are formed non-quasistatically and therefore (a) they 
are launched directly into the non-linear regime of gravitational 
collapse, and (b) they are likely to have retained some internal 
turbulence. 

Burkert & Bodenheimer (2000) have pointed out that the 
internal turbulence in molecular cores can be represented by a 
Gaussian random velocity field having a power spectrum of the 

form P(k) oc k n , with r\ 3 to - 4. This not only reproduces 

the observed scaling relation between size and linewidth. It also 
reproduces the observed rotation of molecular cores. Thus there 
is no need to invoke ordered rotation as an additional source of 
support in molecular cores, and indeed there is no observational 
evidence for significant ordered rotation in prestellar cores (e.g. 
Jessop & Ward-Thompson 2001). 

We have therefore undertaken a numerical study of the 
influence of turbulence on the collapse and fragmentation of 
prestellar cores. We have taken as our reference point a simple 
model of the core L1544, and in the first paper of this series 
(Goodwin et al. 2004a, hereafter Paper I) we have shown that 
cores with even a low level of turbulent energy routinely spawn 
multiple stellar systems. Specifically, in an ensemble of 20 sim- 
ulations of the collapse of 5.4M cores having an initial ratio 
of turbulent to gravitational energy 



ttturb = 



t/turb 

mi 



: 0.05 , 



(1) 



80% of the cores form at least two, and in one case ten, ob- 
jects (stars and brown dwarfs). In addition, the distributions of 
semi-major axis, mass ratio and eccentricity for the resulting 
multiple systems are consistent with the distributions for ob- 
served binary systems (e.g. Duquennoy & Mayor 1991, here- 
after DM91). Paper I has also shown that while low levels of 
turbulence can easily produce multiple fragmentation, the out- 



come of any one simulation depends sensitively on the exact 
details of the turbulent velocity field. Consequently a statistical 
approach is essential in evaluating the influence of turbulence 
on multiple star formation in cores. 

In this paper we extend the simulations of Paper I to ex- 
amine the effect of different levels of turbulence on star forma- 
tion within dense molecular cores. Using exactly the same core 
structure as in Paper I we simulate ensembles of between 10 
and 20 cores with a turb = 0, 0.01, 0.025, 0.05, 0.10 and 0.25. 
We examine the numbers and masses of stars and brown dwarfs 
that form and the properties of the multiple systems to which 
some of them belong. 

We note that these levels of turbulence involve much lower 
non-thermal velocities than the earlier work of Whitworth 
et al. (1995), Turner et al. (1995), Whitworth et al. (1996), 
Klein et al. (2001, 2003), Bate et al. (2002a,b, 2003), Bonnell 
et al. (2003), Delgado-Donate (2003, 2004) Goodwin et al. 
(2004a,b) Hennebelle et al. , and Hennebelle et al. (2003, 
2004). Consequently they may be applicable to scenarios in 
which instability develops more quasistatically due to ambipo- 
lar diffusion, provided that some turbulence can persist through 
(or be regenerated after) the ambipolar diffusion phase, and 
provided the subsequent collapse is sufficiently rapid. 

In Section 2 we define the initial conditions for the simu- 
lations. In Section 3 we describe the code and the constitutive 
physics used. In Section 4 we outline the different ensembles 
of simulations performed with different values of a tur b, and in 
Section 5 we collate the statistics from the different ensembles. 
Section 6 discusses the statistics in terms of the underlying 
physics, and Section 7 gives our main conclusions. 

2. Initial conditions 

Molecular cores which are associated with IRAS sources are 
presumed to have already formed a protostar, and are classified 
as protostellar cores, whereas those which have no associated 
IRAS source are classified as starless cores (Beichman et al. 
1986). The densest starless cores are presumed to be on their 
way to forming stars, and are therefore classified as prestellar 
cores (Ward-Thompson et al. 1994). We base our initial condi- 
tions on the observed properties of prestellar cores. 

The density in a pre-stellar core is approximately uniform 
in the inner few thousand au, but further out it decreases as 
r~ n with 2 ^ 77 $ 5, and eventually it merges with the back- 
ground (e.g. Ward-Thompson et al. 1994; Andre et al. 1996; 
Ward-Thompson et al. 1999; Andre et al. 2000; Tafalla et al. 
2004). A good fit to the density in a pre-stellar core is given by 
a Plummer-like profile, 



Pkernel 



/^kernel) 2 ) 2 ' 



(2) 



where pkemei is the central density and kernel is the radius of 
the region in which the density is approximately uniform. We 
set pkemei = 3 x 10~ 18 gcm~ 3 and /?kernei = 5000 au, with the 
outer boundary of the core at R C0K = 50000 au, so the core has 
total mass M core = 5.4M©. The core is initially isothermal with 
T = 10 K, and uniformly molecular with mean gas-particle 
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mass m = 4 x 10 24 g. This means that the core has a ratio 
of thermal to gravitational energy of 



«the 



tAherm 



0.3 



(3) 



Thus far, these are exactly the same initial conditions as we 
used in Paper I. 

Molecular cores show a significant non-thermal contribu- 
tion to their line-widths, which is attributable to internal turbu- 
lence. Fig.^shows the estimated ratios of turbulent to gravita- 
tional energy, a tur b (see Eqn. {TJ), and the estimated masses, 
M C0K , for prestellar cores from the Jijina et al. (1999) cata- 
logue. These cores were selected as prestellar on the basis of 
their having low temperature (< 20 K), no associated IRAS 
source and no observed outflow. The shaded area in Fig. \l\ 
shows the region of parameter space that the simulations in this 
paper cover, i.e. a 5.4M core with a range of a t urb from to 
0.25. 

To model the turbulence, a divergence-free gaussian ran- 
dom velocity field is superimposed on the core (cf. Bate et 
al. 2002a,b; Bate et al. 2003; Fisher 2004; Bonnell et al. 
2003; Delgado-Donate et al. 2003, 2004). The power spec- 
trum of the velocity field is P(k) cc k~ 4 ', so as to mimic the 
observed scaling between size and line-width in interstellar 
gas clouds (Larson 1981; Burkert & Bodenheimer 2000). The 
magnitude of the velocity field is normalized so that a tur b = 
0, 0.01, 0.025, 0.05, 0.1 or 0.25 , and at least ten realizations 
have been made with each value of a tur b (as summarised in 
Tabledi. 

In each realization the random number seed for the turbu- 
lence is different, and hence the detailed structure of the veloc- 
ity field is different. It is essential that many different realiza- 
tions be performed for a given value of a t urb, because the mix of 
protostars and multiple systems which forms depends critically 
on the details of the velocity field. Therefore different values 
of ff tur b can only be compared in a statistical sense by perform- 
ing an ensemble of different realizations for each representative 
value of aturb- 

Each simulation is evolved for 0.30 Myr. We have chosen 
this end-time for three reasons, (i) It facilitates comparison 
with the a tu ,b = 0.05 simulations reported in Paper I. (ii) In 
most simulations, object formation (i.e. sink creation) ceases 
around 0.15 to 0.20 Myr, and in only one case is a new object 
formed after 0.25 Myr. (iii) By this stage, ~ 66 % of the mass 
is already in stars and brown dwarves (this fraction decreases 
slightly with increasing ffturbX an d so the remaining diffuse gas 
is likely to be affected by negative feedback from the existing 
stars and brown dwarves; feedback is not included in these sim- 
ulations, although we are currently exploring its effect (Boyd et 
al. in preparation). 

The properties of the resulting protostars in each ensemble 
are compared, both as a function of aturb; an d against the ob- 
served properties of young stellar objects. In this latter regard, 
the observational samples used for comparison are the local 
field G-dwarfs studied by DM91 and the multiple systems in 
Taurus studied by White & Ghez (2001; hereafter WG01). 



3. Computational method and constitutive physics 

The simulations are performed using SPH (Lucy 1977; Gingold 
& Monaghan 1977; Monaghan 1992). Our SPH code (dragon) 
uses the standard M4 kernel (Monaghan & Lattanzio 1985) and 
invokes variable smoothing lengths (so that each particle has 
Afneib = 50+5 neighbours). An octal tree is built to facilitate the 
computation of gravitational accelerations and the identifica- 
tion of neighbours. Gravity is kernel-softened with the particle 
smoothing length, and standard artificial viscosity (Monaghan 
1992) is included, with a v = 1 and /3 V = 2. 
We use a barotropic equation of state, 



Pip) _ 2, , 2 



1 + 



p 

Pcrit 



(4) 



Thus at low densities (p < p crrt = 10~ 13 gcirT 3 , where radiative 
cooling is optically thin) the gas is approximately isothermal 
with c s - Co = 0.2 km s -1 , corresponding to molecular gas at 
temperature T = 10 K; and at high densities, (p > p cl -i t , where 
cooling radiation is trapped by the high optical-depth) the gas 
behaves adiabatically, with ratio of specific heats y 5/3. 
(y - 5/3 is the appropriate ratio of specific heats because the 
rotational degrees of freedom of H2 are not excited until the 
temperature rises above ~ 300 K, and our simulations do not 
follow the gas dynamics to such high temperatures; see below). 
Eqn. (0} reproduces the thermal behaviour of the gas in detailed 
one-dimensional simulations of the collapse of 1M cores (e.g. 
Larson 1969; Tohline 1982; Masunaga & Inutsuka 2000). 

Whenever a gravitationally bound condensation forms and 
the density of an SPH particle within the condensation rises 
above 100p cl -i t = 10~" gem -3 , all the particles within 5 au of 
that particle are replaced with a sink particle having radius 
^sink = 5 au. (From Eqn. @, we can estimate that the tempera- 
ture of gas at p = lOOpcrit is ~ 225 K, and so the introduction of 
sink particles makes it unnecessary to treat the thermal behavi- 
ous of the gas at temperatures above ~ 225 K.) Sink particles 
interact with the gas both gravitationally, and by accreting SPH 
particles that enter the sink radius and are bound to the sink 
(see Bate et al. 1995 for a detailed description of sink particles). 
As in Paper I, we refer to sink particles as 'objects', and then 
more specifically to 'stars' when the sink mass is > O.O8M 
and 'brown dwarfs' when the mass is lower than this. 

4. Ensembles of simulations with different <r uirb 

4.1. a tur b - (no turbulence) 

As reported in Paper I, when a core has no turbulence only one 
object forms, very close to the centre of the core. The evolution 
follows very closely the semi-analytic model of Whitworth & 
Ward-Thompson (2001). In particular, the accretion rate is very 
large early on, and then decreases. After 0.3 Myr the stellar 
mass reaches 3.75M . (In this case the ten different realizations 
involve different initial SPH particle positions.) 

4.2. ^urb = 0.01 

The ensemble of ten simulations with aturb = 0.01 is virtually 
identical to that with Cturb = 0.00 (no turbulence). Only one 
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Fig. 1. The filled circles give estimated values of a m ±, ( me rat i° of turbulent to gravitational energy, see Eqn. Q) and M core (core 
mass) for the starless cores in the Jijina et al. (1999) catalogue. The hatched region gives the range of (ofturt» M core ) values used in 
this paper (i.e. < aturb < 0.25 and M core = 5.4M ). The cross (x) and the star (*) give - respectively - the (aturb, M cole ) values 
used by Bate et al. (2002a,b;2003) and Delgado Donate et al. (2003a,b). 



star ever forms, and this happens about one free-fall time (~ 
0.055 Myr) after the start of the simulation, always close to the 
centre of mass of the core. This level of turbulence is apparently 
too low to induce multiple fragmentation. 

4.3. a t urb = 0.025 

Cmrb = 0.025 appears to be approximately the minimum level 
of turbulence required for multiple object formation. Of the ten 
simulations with a t urb = 0.025, eight produce only a single 
star (as with lower levels of turbulence), but one simulation 
produces six objects, and one produces nine objects. 

Specifically, this last simulation produces three 
intermediate-mass stars in an hierarchical triple system 
embedded in the core (a tight binary with component masses 
1.07M Q and 0.88M Q and semi-major axis 8.6 au, plus a third 
star with mass 1.65M Q orbiting at ~ 70 au). In addition, the 
simulation produces one very low-mass star (0.087M G ) and 
five brown dwarfs (O.O34M to 0.072M G ), all of which are 
ejected from the core. 

4.4. a turb = 0.05 

Of the twenty simulations performed with or tur b = 0.05, four 
produce just a single star, and the remaining sixteen produce 
71 stars and 16 brown dwarfs between them (between 2 and 
10 objects per simulation). Of these 71 stars, 44 remain in the 
core in multiple systems, and the rest are ejected from the core. 
Of the 16 brown dwarfs, 15 are ejected, and only one remains 
in a binary system in the core. The mean number of objects 
formed per simulation is 4.55. Further details of this ensemble 
of simulations are given in Paper I. 



4.5. a twb = 0.10 

Of the twenty simulations performed with a tur b = 0.10, five 
produce just a single star, and the remaining fifteen produce 76 
stars and 19 brown dwarfs between them (between 3 and 10 
objects per simulation). Of these 76 stars, 53 remain in the core 
in multiple systems, and the rest are ejected from the core. Of 
the 19 brown dwarfs, 17 are ejected and only 2 remain in mul- 
tiple systems in the core. The mean number of objects formed 
per simulation is 4.75. 

4.6. aturb = 0.25 

In the ten simulations with fftui-b = 0.25, a total of 60 objects are 
produced, 49 stars and 1 1 brown dwarfs, with each simulation 
producing between 3 and 10 objects. Of the 49 stars, 36 remain 
in the core in multiple systems, and the rest are ejected from 
the core. All 1 1 brown dwarfs are ejected from the core. 

5. Statistics 

Details of the numbers, masses and multiplicities of the ob- 
jects produced in each of the fifty simulations with a m ±, = 
0.05, 0.10 and 0.25 are shown in Table |2 The simulations 
with low turbulence (a t ui-b ^ 0.025) are omitted because al- 
most all of them produce only a single object, and therefore 
the discussion will now concentrate on the simulations with 
a tmb =0.05, 0.10 and 0.25. 

5.1. Numbers of objects and formation timescale 

Increasing the level of turbulence has several effects. In the 
first instance, it delays somewhat the time at which objects are 
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Table 1. For each value of a t urb, we list the number of realizations simulated (MeaiX the mean mass in objects at the end of 
the simulations (M tot /M©), the average number of objects formed per simulation ^,/Vobj)), the companion probability (cp), the 
companion frequency (cf), the multiplicity frequency (mf) and the pairing factor (pf). The companion frequency, cf, is given for 
all objects and then seperately for low-mass objects (M < 0.5M G ) and for high-mass objects (M > O.5M ). 



& turb 


Meal 


M tot /M e 


(Km) 


cp 


all 


cf 

low 


high 


mf 


pf 


0.00 


10 


3.75 


i 














0.01 


10 


3.75 


l 














0.025 


10 


3.70 


2.3 ±0.5 














0.05 


20 


3.57 


4.6 ± 0.5 


0.51 


0.95 


0.30 


1.45 


0.27 


1.71 


0.10 


20 


3.35 


4.8 ±0.5 


0.59 


1.39 


0.48 


2.32 


0.32 


2.11 


0.25 


10 


3.14 


6.0 ± 0.5 


0.62 


1.57 


0.81 


2.38 


0.32 


2.36 




0.01 



Fig. 2. The average number of objects formed in a core, ^A/obj), as a function of the ratio of turbulent to gravitational energy orturb- 
The error bars give the standard deviation (see Eqn. ([5}. 



formed, and at the same time it increases the average number 
of objects formed. 

With ff tu ,b = 0.05, the first object (hereafter the primary 
protostar) forms 0.05 Myr to 0.06 Myr after the start of the 
simulation, and ~ 95% of all the other objects have formed 
by 0.12 Myr. With or tur b = 0.25, the primary protostar forms 
0.06 Myr to 0.08 Myr after the start of the simulation, and ob- 
jects continue forming up to 0.15 Myr; this is because the extra 
turbulent energy gives the core extra support, and therefore de- 
lays its collapse. 

The majority of secondary objects form in a dense disc- 
like slab around the primary protostar. The instabilities which 
produce these secondary objects are usually seeded - and pro- 
pelled into the non-linear condensation regime - by the inho- 
mogeneities in the accretion flow onto the slab. Higher levels of 



turbulence result in larger inhomogeneities, and hence in more 
secondary objects. 

For a tur b < 0.01 the mean number of objects formed is 
(Afobj) = 1 (i-e. there are no secondary objects), whereas for 
ffturb = 0.25, (N bj) = 6 (i.e. there are on average five sec- 
ondary objects). Values of or tur b and (A/'obj) are listed in Table 
1, and plotted on Fig. f2] where the error bars give the standard 
deviation, i.e. 



(O-^obj) 2 



1/2 



(5) 
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Table 2. A summary of the results of the simulations with a tur b = 0.05, 0.10 and 0.25 , at time t = 0.3 Myr. Column 1 gives the 
simulation identifier and Column 2 gives a tur b. Column 3 gives M Q bj, the total mass of objects formed (stars plus brown dwarfs), 
Column 4 gives W bj> the total number of objects formed, and Column 5 gives A/bd> the total number of brown dwarfs formed. 
Column 6 gives the multiplicities of the multiple systems formed, and the final column (Column 7) gives the mass of each 
individual object. Those objects which are part of a binary system are distinguished with b , those which are part of a triple system 
are distinguished with a ', and those which are part of a quadruple system (or in one case a quintuple system) are distinguished 
with q . 



ID 


ftturb 


M ob] 


A'obj 




Multiplicity 


Masses/M 


071 


0.05 


2.94 


7 


2 


Triple 


1.31', 0.61', 0.52', 0.27, 0.12, 0.063, 0.048 


072 


0.05 


3.72 


4 





Binary 


2.32*, 0.74*, 0.48, 0.18 


073 


0.05 


3.10 


10 


2 


Binary + Triple? 


1.07", 0.66", 0.43, 0.34, 0.17, 0.13', 0.10', 0.09', 0.076, 0.040 


074 


0.05 


4.02 


3 





Triple 


1.63', 1.56', 0.83' 


075 


0.05 


3.69 


2 





Binary 


2.63*, 1.06* 


076 


0.05 


3.61 


3 


1 


Binary 


2.18*, 1.40*, 0.028 


077 


0.05 


3.75 


6 


1 


Triple 


1.60', 1.16', 0.64', 0.18, 0.12, 0.050 


078 


0.05 


3.65 


7 


2 


Triple 


1.09', 1.03', 0.69', 0.58, 0.18, 0.045, 0.041 


079 


0.05 


3.81 


8 


3 


Triple 


1.27', 1.16', 0.69', 0.39, 0.21, 0.044, 0.030, 0.025 


080 


0.05 


3.63 


1 





Single 


3.63 


081 


0.05 


3.69 


1 





Single 


3.69 


082 


0.05 


4.01 


4 





Quadruple 


1.52 s , 0.91", 0.89 ? , 0.69 9 


083 


0.05 


3.56 


4 





Triple 


1.43', 0.83', 0.70', 0.60 


084 


0.05 


3.55 


5 





Binary 


1.46*, 1.28*, 0.43, 0.19, 0.18 


085 


0.05 


3.47 


8 


3 


Binary 


1.43*, 0.76°, 0.51, 0.47, 0.14, 0.064, 0.045, 0.039 


086 


0.05 


3.94 


7 


1 


Triple 


1.23', 1.03', 0.73, 0.71', 0.11, 0.098, 0.027 


087 


0.05 


3.67 


2 





Binary 


3.19*, 0.48* 


088 


0.05 


3.35 


1 





Single 


3.35 


089 


0.05 


3.61 


7 


1 


Quadruple 


1.20* , 0.89 9 , 0.57, 0.51, 0.29 ? , 0.11, 0.041'' 


090 


0.05 


2.62 


1 





Single 


2.62 


001 


0.10 


3.78 


3 





Triple 


1.49 , 1.15 , 1.13 


002 


0.10 


2.83 


1 





Single 


2.38 


003 


0.10 


3.72 


1 





Single 


3.72 


004 


0.10 


3.48 


1 





Single 


3.48 


005 


0.10 


2.86 


4 


1 


Binary 


1.43*, 0.65*, 0.77, 0.02 


006 


0.10 


2.84 


1 





Single 


2.84 


007 


0.10 


3.15 


5 





Triple & Binary 


1.76', 0.72', 0.47', 0.10*, 0.10* 


008 


0.10 


3.22 


6 


2 


Quadruple 


1.97 s , 0.47 ? , 0.35'', 0.34 9 , 0.03, 0.06 


009 


0.10 


3.48 


8 


4 


Quadruple 


2.28' ? , 0.49 ? , 0.26 ? , 0.25'', 0.05, 0.08, 0.04, 0.04 


010 


0.10 


3.31 


8 


1 


Quadruple 


0.76 9 , 0.74 9 , 0.58 ? , 0.57' 7 , 0.08, 0.09, 0.03, 0.46 


011 


0.10 


3.96 


12 


4 


Triple & binary? 


0.89', 0.82', 0.82', 0.04*, 0.04*, 0.42, 0.38, 0.03, 0.03, 0.25, 0.12, 0.11 


012 


0.10 


3.60 


6 


2 


Triple 


1.34', 0.92', 0.79', 0.50, 0.04, 0.02 


013 


0.10 


3.18 


10 


3 


Quadruple & binary 


0.77'', Q.6% q , 0.61 9 , 0.60?, 0.11*, 0.11*, 0.10, 0.05, 0.04, 0.06 


014 


0.10 


3.29 


4 


1 


Binary 


1.58 , 1.16 , 0.49, 0.08 


015 


0.10 


2.48 


1 





Single 


2.48 


016 


0.10 


3.58 


4 





Triple 


1.23 , 1.15 , 1.11', 0.09 


017 


0.10 


3.41 


8 





Quintuple 


1.10 s , 0.98' 7 , 0.32' 7 , 0.27 ? , 0.27'', 0.15, 0.14, 0.17 




n 1 n 


? A Q 




A 

u 


Quadruple 




019 


0.10 


3.58 


5 


1 


Triple 


1.38', 1.03', 1.00', 0.11,0.06 


020 


0.10 


3.77 


3 





Triple 


1.28', 1.27', 1.22' 


041 


0.25 


3.13 


6 





Quadruple 


0.79«, 0.71", 0.68", 0.36", 0.33, 0.27 


042 


0.25 


3.14 


8 


2 


Triple & binary 


1.92', 0.41', 0.30', 0.10*, 0.18*, 0.07, 0.14, 0.02 


043 


0.25 


2.69 


5 


1 


Triple 


0.81', 0.81', 0.53', 0.50, 0.03 


044 


0.25 


3.17 


5 





Quadruple 


1.26", 0.92", 0.46", 0.27", 0.26 


045 


0.25 


3.18 


8 


2 


Quadruple 


0.78", 0.55", 0.54", 0.47", 0.05, 0.04, 0.62, 0.13 


046 


0.25 


3.19 


3 





Binary 


1.70*, 0.94*, 0.54 


047 


0.25 


3.31 


5 





Quadruple 


1.07", 0.81", 0.64", 0.46", 0.33 


048 


0.25 


3.36 


10 


4 


Quadruple 


0.87", 0.76", 0.56", 0.56", 0.08, 0.03, 0.02, 0.06, 0.37, 0.10 


049 


0.25 


2.84 


3 





Triple 


2.01', 0.42', 0.41' 


050 


0.25 


3.37 


7 


2 


Quadruple 


0.91", 0.91", 0.68", 0.67", 0.02, 0.06, 0.12 
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Fig. 3. The normalized mass functions for a lm \, = 0.05 (top), 0. 10 (middle) and 0.25 (bottom). The shaded regions show stars 
stable multiple systems, the hashed regions show stars in unstable multiple systems, and the open regions show single stars. 
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Fig. 4. The combined (and un-normalized) mass function for all the simulations having a- tul -b > 0.05, (i.e. aturb = 
0.05 (x20), 0.10 (x20) and 0.25 (x 10)), this being an approximation to the distribution of ff tur b values in the Taurus star formation 
region. The error bars give + uncertainties. The dotted lines show separate fits to the unbound objects (the flat distribution at 
low mass) and the bound objects (the log-normal distribution at high mass); parameter fits are given in the main text). The dashed 
line shows the combination of these two separate fits. 



5.2. The mass function 

The normalized mass functions (MF) for the ensembles with 
aturb = 0.05, 0.10 and 0.25 are shown in Fig. [3] In each case 
the shaded region shows the MF of objects in multiple systems 
and the open region shows the MF of single objects. The hashed 
region in the top figure shows an unstable triple that formed 
late on in one simulation, and is therefore likely to decay into a 
binary and an ejected single (see Paper I for more details). 

The MFs are clearly very similar for each value of arturb> viz. 
a high-mass peak of predominantly multiple stars and a low- 
mass tail of ejected stars and brown dwarfs. However, as a tur b 
increases, the MF shifts slightly to lower masses. There are two 
reasons for this, (a) For higher aturb the overall collapse is de- 
layed by the extra turbulent support, and therefore when the 
simulations are terminated at 0.3 Myr less mass has been incor- 
porated into objects (the mean mass incorporated into objects 
is given in the third column of Table 1). (b) For higher a tul b, 
the accretion flow onto the disc around the primary protostar is 
lumpier, so more objects are formed but individually they are 
less massive. 

The combined mass function for all simulations with 
«turb ^ 0.05 is shown in Fig.0] The distribution of high-mass 
(predominantly bound) stars is well fitted by a log-normal dis- 
tribution having mean {£ogu)[M}) = 0.05 and standard devia- 
tion <X( ogm [M] = 0.04. The distribution of low-mass (predom- 
inantly unbound) objects is consistent with being flat in log- 
space from our resolution limit at ~ O.O25M up to ~ O.5M , 
above which it declines. The high-mass (predominantly bound) 
stars have an average mass of ~ 1M because, after ejections 
have removed some objects, there are usually two to four stars 
left in the core, and they are then able to accrete a total of 



~ 3M between them. A more massive core would spawn more 
massive stars (Goodwin et al. in preparation). 

The fraction of objects which are brown dwarfs is 
A/BD/Afobj ~ 18%, and there does not appear to be a systematic 
dependence on the level of turbulence. This is somewhat higher 
than in Taurus (Nbd/ 'Afobj ~ 13%; Briceno et al. 2002), and 
somewhat lower than in Orion (AfBD/Af bj ~ 26%; Muench et 
al. 2002). The fraction of low-mass objects (M < O.5M ) also 
appears to be independent of a tm b, and approximately ~ 50%. 

5.3. Companion star frequencies 

In analyzing the multiplicity of the objects formed in our simu- 
lations, we define 'systems' to include single objects, and 'mul- 
tiple systems' to include only systems containing more than 
one object. The primary is the most massive star in a system; 
in a single it is the only star. Thus, if S is the number of single 
objects and B, T, Q and Q' are the numbers of binary, triple, 
quadruple and quintuple systems, respectively, the total num- 
ber of objects is (S + 2B + 3T + 4Q + 5Q + ...), the total number 
of systems is (S + B + T + Q + Q' + ...) (which is the same as 
the total number of primaries), and the total number of multiple 
systems is (B + T + Q + Q + ...). 

Many different statistics have been introduced as mea- 
sures of stellar (and brown dwarf) multiplicity (e.g. Reipurth & 
Zinnecker 1993), and they all reflect slightly different things 3 . 
They can be divided into two groups. 



3 Unfortunately, the different measures usually have more than one 
name. For consistency we have adopted the nomenclature proposed by 
Reipurth & Zinnecker 1993. 
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The first group of measures is normalized to the total num- 
ber of objects, and is useful because it is straightforward to 
derive these measures for a subset of objects (for example, low- 
mass stars, in the range 0.08M Q < M < O.5M ). The compan- 
ion probability (Reipurth & Zinnecker 1993), 



cp 



2B + 3T + 4Q + 5Q' + ... 
S + 2B + 3T + 4Q + 5Q' + . 



(6) 



gives the probability of an object having at least one compan- 
ion, or equivalently the fraction of objects which is in multi- 
ple systems. However, it gives no indication of whether objects 
with companions are in binaries, or triples, or higher multiples. 
Therefore we prefer the companion frequency, 



cf = 



2B + 6T + 12(2 + 202' + • 
S +2B + 3T + 4Q + 5Q' + 



(7) 



which gives the mean number of companions per object. 

The second group of measures is normalized to the total 
number of systems or the total number of multiple systems. 
The multiplicity frequency (RZ93), 



mf 



B + T + Q + Q' + ... 
S + B + T + Q + Q + ... 



(8) 



gives the fraction of systems which is multiple. The pairing 
factor (RZ93), 



Pf 



B + 2T + 3Q + AQ' + 
B + T + Q + Q' + ... 



(9) 



gives the mean number of orbits per multiple system, or equiva- 
lently the mean number of companions per primary in multiple 
systems. 

Table ^ records the values of all these measures for cores 
having a tur b = 0.05, 0.10 and 0.25. We have also calculated the 
companion frequency, cf, separately for the low-mass objects 
(< O.5M ) and the high-mass objects (> O.5M ). The fraction 
of objects in multiple systems (i.e. the companion probability, 
cp) increases steadily with increasing a m \,. The mean num- 
ber of companions (i.e. the companion frequency, cf) increases 
even more rapidly with increasing ar tur b. cf is always larger 
for high-mass objects than for low-mass objects (i.e. there are 
fewer low-mass objects in multiples than high-mass objects), 
but the rate of increase of cf with increasing a tur b is greater for 
the low-mass objects. 

The smaller companion frequency for low-mass objects is 
due to the fact that low-mass objects have usually been ejected 
from the core before they could accrete much mass (that is why 
they have low mass), and ejected stars tend to be singles. Stars 
that end up in stable multiples also tend to remain in the core, 
and therefore they grow to larger masses by continuing to ac- 
crete. 

The multiplicity of low-mass stars increases with increas- 
ing arturb. because a greater number of higher-order multi- 
ples is formed in cores with higher o-tmb- For example, when 
ffturb = 0.05, only 10% of simulations produce quadruples, but 
this fraction rises to 60% for a- tur b = 0.25. In higher-order mul- 
tiples, the low-mass objects tend to be outlying members. They 
are less able to accrete from the remaining gas, and they tend 
to remain low-mass. 



5.4. Separations 

Fig. |3 shows the cumulative distribution functions (CDFs) of 
the semi-major axes of multiple systems, for different initial 
levels of turbulence, a tur b = 0.05, 0.10 and 0.25. For com- 
parison the gaussian fit to the DM91 local G-dwarf sample is 
plotted as a dashed line. 

As noted in Paper I, the semi-major axis distribution for 
the of tur b = 0.05 ensemble is consistent with the DM91 obser- 
vations. In contrast, the semi-major axis distributions for the 
ffturb = 0.10 and 0.25 ensembles both have too many hard bi- 
naries (a < 20 au), and both are rejected by the KS test as be- 
ing drawn from the DM91 fit, at > 90% confidence. A similar 
excess of hard binaries is predicted by the core fragmentation 
simulations of Delgado-Donate et al. (2003, 2004), which in- 
voke even higher levels of turbulence (arturb = 1), and by the 
iV-body simulations of Sterzik & Durisen (2003). 

As described in Paper I, hard binaries are formed primarily 
by few-body interactions, including those which eject low mass 
objects. Consequently, in simulations where larger numbers 
of objects are formed, the binaries are on average harder (re- 
gardless of a tur b). For example, when few objects are formed, 
say A^,bj < 3, the average separation of binaries is > 100 au, 
whereas when N bj = 4 the average separation is ~ 30 au, and 
when iV bj > 5 it is ~ 20 au. Fig.[6]shows the semi-major axes 
of all systems plotted against the number of objects formed in 
that simulation. There is clearly a trend of decreasing semi- 
major axis with increasing number of objects. 

Since the gravitational forces between objects are kernel 
softened with a smothing length h equal to the sink radius 
^sink = 5 au, orbits with small semi-major axes (a < 5 au) are 
also softened. It follows that the distribution of semi-major axes 
below 5 au is distorted. Given that the code conserves angular 
momentum very accurately, we infer that these already hard or- 
bits should be even harder. For a tm \, =0.10 and 0.25 this would 
exacerbate the difference between the numerically derived dis- 
tribution of semi-major axes and the observations of DM91. 
Conversely, for arturb - 0.05, it would improve the agreement 
with the DM91 distribution. 

5.5. Mass ratios 

Fig.0shows the distribution of binary mass ratios, q = MijM\, 
for all simulations having a tm \, > 0.05. Note that high q means 
q ~ 1, i.e. components of comparable mass. 

For of tur b = 0.05, the distribution of mass ratios is quite flat, 
and reminiscent of the observed distribution for local G dwarfs 
(DM91;Mazehet al. 1992). 

For a = 0.10 and 0.25, the distribution is dominated by 
high-g close binaries. All binaries in these ensembles have 
semi-major axes a < 100 au (the high-a tail is produced by 
wider orbits in higher-order systems), and 64% of these have 
q > 0.8. This is very similar to the mass ratio distribution ob- 
served in Taurus- Auriga by WG01, who found that over ~ 60% 
of binaries with separations < 100 au had q > 0.8. In our sim- 
ulations, systems with high mass ratio tend to be close (all 
systems with q > 0.7, and most systems with q > 0.4, are 
binaries with a < 20 au), but the reverse is not always true: 
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Fig. 5. The histograms show the cumulative distribution functions of semi-major axes for the ensembles with a tur b = 
0.05 (top), 0.10 (middle) and 0.25 (bottom). The dashed line shows the guassian fit to the DM91 period distribution. 
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Fig. 6. For each binary system, the semi-major axis, a is plotted against the number of objects formed in that simulation, N bj, 
with open circles for for o- tur b = 0.05, filled circles for a t urb = 0.10, and stars for a tul -b = 0.25. The average semi-major axis, 
a, decreases with increasing N \,j, due to dynamical hardening. To the left of the dashed line at a = 5 au the a-values are upper 
limits, due to gravity softening. 
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Fig. 7. The distribution of mass ratios for simulations with a tul -b = 0.05 (filled), a- tur b = 0.10 (hashed) and a tur \, = 0.25 (open). 
(Note that the number of simulations in the a tur b = 0.25 ensemble is only half of the number in the other ensembles.) 



in other words, there are a few close binaries with low mass 5.6. Eccentricities 
ratios. Close binaries are presumed to acquire high mass ra- 
tios because the material accreting onto the system has high 
specific angular momentum, and is therefore more readily ac- 
commodated by the secondary (Whitworth et al. 1995; Bate & 
Bonnell 1997; Paper I). 



Fig. [8] shows the CDF of eccentricity for all the simulations 
having a tul -b > 0.05 and the linear fit to the observed distribu- 
tion proposed by DM91. The two distributions are consistent. 
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Fig. 8. The cumulative distribution function of eccentricities for all simulations with a tur b > 0.05. The dashed line shows the 
DM91 fit of a linearly increasing eccentricity distribution. 



6. Discussion 

6.1. Minimum level of turbulence for multiple formation 

When the initial level of turbulence is low, a lm \, < 0.01, it 
seems that a core can only spawn a single central star. Even 
for ff tu ,b = 0.025, the core is unlikely to spawn a multiple sys- 
tem. We therefore focus our discussion on the higher levels of 
turbulence, a t urb = 0.05, 0.10 and 0.25, for which multiple star 
formation is the norm. In this range a number of significant 
systematic trends are evident. 

6.2. Time-scale for star formation 

As Cturb increases from 0.10 to 0.25, the average timescale for 
star formation increases somewhat, due to the extra support 
which turbulence affords the core. For arturb = 0.05, the pri- 
mary protostar forms after 0.05 Myr, and most of the secondary 
protostars have formed by 0.12 Myr. For ar tur b = 0.10 and 0.25, 
the primary protostar forms after 0.06 Myr, and most of the sec- 
ondary protostars have formed by 0.15 Myr. In only one case 
(run 073) do objects form after 0.25 Myr, and so it appears that 
the fragmentation phase is almost always over by the end of the 
simulations at 0.3 Myr. After the end of the simulations accre- 
tion will be on-going. However, feedback from the protostars is 
expected to become very important, possibly dispersing a sig- 
nificant fraction of the gas not already in stars or discs around 
them. 

6.3. Number of objects formed 

As a tur b increases from 0.10 to 0.25, the average number of 
objects formed, ^W bj), increases from 4.6 + 2.6 to 6.0 ± 2.0 
(see Fig. This is because a higher level of turbulence gen- 



erates more density contrast - i.e. more numerous and more 
compressed lumps - and therefore more protostars. 

6.4. Masses of objects formed 

As a t urb increases from 0.10 to 0.25, the average mass of the 
objects formed decreases slightly (see Fig. [3J|. There are two 
factors involved here. First, as noted above a higher level of tur- 
bulence means that more objects are formed. Second, a higher 
level of turbulence means that the core has more support, and 
therefore a smaller fraction of its mass has condensed out after 
0.3 Myr. 

6.5. The bimodal mass function 

Apart from this slight decrease in average mass with increasing 
turbulence, the form of the mass function appears to be inde- 
pendent of arturb. Specifically, the mass function is bimodal: the 
lower-mass stars (which tend to be single stars ejected from the 
core) subscribe to a flat segment of the mass function; and the 
higher-mass stars (which tend to be those remaining in the core 
and pairing up in multiple systems) subscribe to a Gaussian 
segment of the mass function (see Fig. @}. 

The critical mass seperating the two modes in the mass 
function, M cl -i t ~ O.5M , arises because of the interplay be- 
tween ejection by dynamical interaction and growth by accre- 
tion. Anosova (1986) has shown that the decay time for small- 
N systems is 

fdecay- 100f cross ~ 17 (— ) l — l yr, (10) 

and the ensemble of low-mass secondaries which forms in 
our cores typically has R ~ 200 au and M ~ 1 to 2M Q , so 
fdecay ~ 3 to 5 x 10 4 yr. The accretion rate is ~ lO _5 M yr _1 , 
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and so the objects which get ejected have masses 5 O.5M . 
The probability of ejection is only weakly dependent on mass, 
oc M~'/ 3 (Anosova, 1986), and so the ejected objects have a flat 
mass function. In contrast, the two or three objects which sur- 
vive the dynamical decay phase remain embedded in the centre 
of the core and compete for the gas which continues to fall into 
the centre, so they grow to ~ M Q . 

6.6. Companion star frequency 

As a tur b increases from 0.10 to 0.25, the companion star fre- 
quency increases slightly for intermediate-mass stars (O.5M 
to 5M ), and quite markedly for low-mass objects (< O.5M ). 

6. 7. Distribution of semi-major axes 

For Qfturb = 0.05, the distribution of semi-major axes is broad 
and indistinguishable from the distribution inferred for local 
G-dwarfs by DM91. There is a lack of very close systems with 
a < 5 au. This is due to the fact that sinks have a finite size 
and their mutual gravity is softened; therefore our code cannot 
resolve very close systems. 

In contrast, for ff tu ,b = 0.10 and 0.25, there are many more 
close systems (5 au ^ a ^ 20 au) than in the DM91 sample (see 
Fig-El, an d this discrepancy would not be alleviated if the code 
were able to resolve very close systems. 

Much of the hardening which produces close binaries is due 
to dynamical interactions with other objects, in particular with 
the low-mass objects which get ejected in the process. The ex- 
cess of close systems produced by higher levels of turbulence 
may therefore be due to the greater number of objects formed, 
and hence the greater potential for dynamical interactions, as 
suggested by Fig. [6] 

6.8. Distribution of mass-ratios 

For arturb = 0.05, the distribution of mass ratios is flat and indis- 
tinguishable from the distribution for local G-dwarfs reported 
by DM91. In contrast, for a t urb = 0.10 and 0.25, there is an 
excess of systems having high mass-ratio, i.e. components of 
comparable mass (see Fig. [7}. Many of these systems with high 
mass-ratio arise because the binary system has accreted mate- 
rial with relatively high specific angular momentum, and this 
material can more easily be accommodated by the secondary 
(e.g. Whitworth et al. 1995). 

6.9. Close systems with comparable components 

For ff tur b = 0.10 and 0.25, the systems with high mass-ratio 
tend also to be close, and it is this sub-population of high- 
mass-ratio close binaries which is the main difference between 
the distributions of semi-major axis and mass-ratio for the pro- 
tostars formed in these simulations, and the distributions of 
semi-major axis and mass-ratio for local G-dwarfs as reported 
by DM91. A similar excess of close systems with comparable 
components was found by Delgado-Donate et al. (2003, 2004), 



who simulated the collapse and fragmentation of cores with 
even higher levels of turbulence (ff tu ,b = 1). 

Taken at face value, this suggests that the local population 
of G-dwarfs must have been formed in cores with low turbu- 
lence (aturb ~ 0.05). However, this conclusion rests on the as- 
sumption that the spherically symmetric 5.4M core and the 
P k oc k~ 4 turbulence spectrum which we have adopted, are rep- 
resentative of the cores forming G-dwarfs, and there is no firm 
basis for this assumption. An alternative explanation is that a 
significant population of close, high-mass-ratio systems has es- 
caped detection, but we believe this to be unlikely. 

A significant contrast to this is found in Taurus. Here WG01 
find that binaries in the separation range lOau S5 a 100 au 
do indeed have significantly higher mass ratios than wider bi- 
naries. Therefore they are compatible with formation in cores 
having higher levels of turbulence, a- tul -b = 0.10 to 0.25. We 
have discussed the origin of the mass function and the binary 
statistics in Taurus in Goodwin et al. (2004b). 

7. Conclusions 

We have explored the influence of turbulence on the fragmenta- 
tion of dense molecular cores, by means of a large ensemble of 
simulations. In this ensemble, we consider a spherically sym- 
metric 5.4M core with a Plummer-like density profile; this is 
a good representation of observed cores like L1544. We seed 
the core with a turbulent velocity field having power spectrum 
P(k) oc k , The number of objects that forms, and the proper- 
ties of the resulting multiple systems depend both on the level 
of turbulence a t uib, and on the details of the turbulent veloc- 
ity field. Therefore for each value of a tur b we we have simu- 
lated many different realizations by changing the random num- 
ber seed for the turbulent velocity field. The main conclusions 
are 

(i) The formation of multiple systems requires a- tul -b <; 

0. 025; a core with a t urb = 0.025 has a ~ 20% chance of forming 
a multiple system, and a core with a- tul -b 0.25 almost always 
forms a multiple system. 

(ii) As a-tm-b is increased, the average time-scale for object 
formation increases, the average number of objects formed in- 
creases, the companion frequency increases (particularly for 
the lower-mass objects), and the average mass of objects de- 
creases. 

(iii) The mass function has a bimodal form. The low- 
mass objects, which are usually single because they have been 
ejected from the core before they could grow above O.5M , 
subscribe to a flat segment of the mass function. The high-mass 
stars, which have usually stayed embedded in the core, grown 
by accretion and paired up with one another in multiple sys- 
tems, subscribe to a Gaussian segment of the mass function. 
Typically 20% of objects are brown dwarfs (M < O.O8M ), 
and 50% are low-mass stars (0.08M o < M < 0.5M Q ). 

(iv) For a tul b ^ 0.10, there is a significant subpopulation 
of binary systems having small semi-major axes and high mass 
ratios, i.e. close systems with components of comparable mass. 
This subpopulation is also found in the simulations of Delgado- 
Donate et al. (2003, 2004) who treat the extreme case a tul -b = 

1 . It is not present in the sample of local G-dwarfs observed 
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by DM91, but there is some evidence for it in the Taurus pre- 
Main-Sequence sample observed by WG01. 

(v) The ensemble of simulations for cores with a tu ,b = 
0.05 reproduces the binary statistics of the DM91 sample 
(companion-star frequency and distributions of semi-major 
axis, eccentricity and mass-ratio) very well. Therefore if the 
other core parameters, which we have not varied (e.g. mass), 
are representative of the cores forming local G dwarfs, we in- 
fer that these cores must have had finite but low levels of tur- 
bulence, O-turb ~ 0.05. 

(vi) Both the mass function, and the binary statistics, for the 
WG01 sample of pre-MS stars in Taurus are reproduced by a 
mix of simulations with a tu ,b = 0.05 (20 realizations), 0.10 (20 
realizations), and 0.25 (10 realizations), as shown by Goodwin 
et al. (2004b). 
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